
Calhoun 

Iniiiiuiiortfl Archive of (he Navjl Pwigndualt School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1990-09 

Aerodynamic coefficients of a symmetrical 
airfoil in an oscillating flow 

Chamberlain, George A. 

Monterey, California: Naval Postgraduate School 


http://hdl.handle.net/10945/10801 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p ://w w w. n ps.e-du/l ib ra ry 


Calhoun is the Naval Postgraduate School's public access digitaI repository for 
research materials and institutional publications created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed —and published —scholarly author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 






NAVAL POSTGRADUATE SCHOOL 
Monterey, California 


r 1 \ 

i 







OCT 0 3 193 

D 




AD-A241 073 


THESIS 


AERODYNAMIC COEFFICIENTS OF A 
SYMMETRICAL AIRFOIL IN AN 
OSCILLATING FLOW 

by 

LT George A. Chamberlain, Jr. 

September, 1990 

Thesis Advisor: James A. Miller 

Approved for public release; distribution is unlimited. 


91-12191 

I fill I lllli lllia .....__ " ■ 








SECURITY ClASS.P'CAT.OiN, Of ThiS PAGE 


REPORT DOCUMENTATION PAGE 

Form Approved 

OMB No 0704 0188 

la REPORT SECURDY C.ASSiF,CATION 

UNCLASSIFIED 

lb RESTRICTIVE MARKINGS 


2b DECLASSIFICATION.'DOWNGRADING SCHEDULE 

A PERFORMING ORGAN.ZAT'ON REPOR' NUMBE 

RiS; 

Naval Postgraduate School 


6a NAME OF PERFORMING ORGAN ZAT.ON 

6 b OFFICE SYMBOL 


(If applicable) 

Naval Postgraduate School 

Code 31 

6c ADDRESS (City, State, and ZIP Code) 


Monterey, CA 93943-5000 



Approved for public release; distribution 
is unlimited. 


S MONITORING ORGANIZATION Rtf OR' *, w V?t ; 


Naval Postgraduate School 


ADD^VCify State and/'PCodr) 
Monterey, CA 93943-5000 



8a NAME Of FUNDING SPONSOR.NG 
ORGANIZATION 


8c ADDRESS (C. ty, State, and ZIP Cod*) 



.Vijf- UN T 
ACCESS 0N N O 


11 T:Nt (Include Security Class firar.uni 

AERODYNAMIC COEFFICIENTS OF A SYMMETRICAL AIRFOIL IN AN OSCILLATING FLOW (UNCLASSI 


12 PERSONA. A . THOR.'?, 

Chamberlain, George A., Jr. 


13a T CPF O' REPORT '3b IMF COVERED 

Master;s Thesis from - 


■A date of REPORT [Year, Month Day) 
September, 1990 



16 SJPP-tVEN'AOY NOTATION 

The views in this thesis are these of the author and do not reflect the official policy 
or position of the Department of Defense or the U. S. Government 


18 SUBJECT TERMS 'Continue on reverse if necessary ana identity by block number) 

Oscillating Flow, Aerodynamic coefficients, Symmetrical 
airfoil 


19 ABSTRACT Continue on reverse if necessa-, and identify by block number) 

A symmetrical airfoil at fixed angles of attack from 0 to 35 degrees was placed in 
an air flow with controlled velocity oscillations at 10 Hz. The instantaneous pressure 
distribution and flow velocity were measured simultaneously. Pressure information was 
integraced numerically to determine instantaneous values for lift, drag and moment about 
the quarter chord. The results were plotted as three-dimensional surface plots in angle of 
attack, cycle time and performance coefficient axes, or with chordwise position, cycle 
time and pressure coefficient axes, as appropriate. The pressure distribution over to 
airfoil was shown to be primarily a function of freestream static pressure. Aerodynamic 
coefficients as a function of angle of attack and time exhibited seemingly random, erratic 
behavior, whereas the mean values of the same coefficients were similar in nature to those 
obtained in steady flow. The data indicate that the instantaneous performance of an 
airfoil in an oscillating flow is not merely a continuum of many quasi-steady states. 


2 0 DS T R E ' ON A v A . A B , ' y O * ABS TRAC¬ 
ED ,\( . AST'S ED D □ SAVE AT «•* (□ o~ r ;SE U S 

2 ' A B j ’ P A r T <,-Cj pTv C \ O'. 

UNCLASSIFIED 

,‘ta NAVE 0 [ 5 iyu'., B.i iNJ 

James A. Miller 

r r O'.; (mclud* At+j Led*-} 

(408) 646-2897 

AAUO 











































Approved for public release; distribution is unlimited. 


Aerodynamic Coefficients of a 
Symmetrical Airfoil in an 
Oscillating Flow 

by 

George A. Chamberlain, Jr. 

Lieutenant, United States Navy 
B.S., United States Naval Academy, 1982 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 

from the 


NAVAL POSTGRADUATE SCHOOL 
September, 1990 



11 









ABSTRACT 


A symmetrical airfoil at fixed angle of attack from 0 to 35 degrees was placed in 
an air flow with controlled velocity oscillations at 10 Hz. The instantaneous pressure 
distribution and flow velocity were measured simultaneously. Pressure information was 
integrated numerically to determine instantaneous values for lift, drag and moment about 
the quarter chord. The results were plotted as three-dimensional surface plots in angle 
of attack, cycle time and performance coefficient axes, or with chordwise position, cycle 
time and pressure coefficient axes, as appropriate. The pressure distribution over the 
airfoil was shown to be primarily a function of freestream static pressure. Aerodynamic 
coefficients as a function of angle of attack and time exhibited seemingly random, erratic 
behavior, whereas the mean values of the same coefficients were similar in nature to those 
obtained in steady flow. The data indicate that the instantaneous performance of an 
airfoil in an oscillating flow is not merely a continuum of many quasi-steady states. 
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I. INTRODUCTION 


The instantaneous pressure distribution on an airfoil in an oscillating flow is not 
well understood at present, chiefly because previous attempts to map it were hampered 
by the inability to measure pressure and velocity simultaneously at data rates compatible 
with the oscillating flow frequency, by the lack of computer resources for data reduction 
and by the lack of graphics software necessary to effectively display the results. The 
widespread use of microcomputers in aerodynamic research and the availability of three- 
dimensional plotting software now makes possible the acquisition, calculation and display 
of the pressure distribution and performance coefficients, C L , C D and C M c/4, in real-time 
format. This t v > pe of information is of particular interest to designers of helicopter rotor 
blades, as it reflects the varying performance of such a wing during its normal cyclic 
operation. 

The basic approach taken in this investigation was to observe the pressure 
distribution over a stationary symmetrical airfoil at angles of attack ranging from 0 to 20 
degrees, in increments of 5 degrees, in a steady flow, and at angles of attack ranging from 
0 to 35 degrees, in increments of 5 degrees, in a controlled oscillating flow at 10 Hz. 
The pressure distribution was converted to coefficient of pressure, C P , and integrated to 
obtain the normal force coefficient, C N . This was further resolved into lift and drag 
components, C L and C D . Additionally, the C P distributions were integrated to obtain the 
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moment coefficient about the quarter chord, C M c/4. This information was then plotted as 
a function of cycle time, t/T, and u.ordwise position, x/c, or as a function of angle of 
attack and cycle time, as appropriate, and a comparison was made of the resulting 
similarities and differences between steady flow and oscillating flow aerodynamic 
coefficients. 





II. EXPERIMENTAL EQUIPMENT 


A. WIND TUNNEL 

1. General Description 

The experiment was conducted in the Naval Postgraduate School Aeronauticai 
Engineering low-speed oscillating flow wind tunnel. Figure 1, a plan view of the wind 
tunnel, shows that it is a continuous flow, in-draft facility. 

The inlet is an eight foot square cross section, which contracts to a two foot 
square cross section in the tes: section. The inlet is equipped with three high solidity 
screens. These, and the 16:1 contraction ratio result in a very smoo’h, steady flow 
through the test section. Measurement of steady flow test section velocity profiles have 
shown that velocity variation is less than 1.0 percent of the mean to within 3 inches of 
any wall. Typical freestream turbulence intensity in steady flow is 0.3 percent. 

The flow is produced by two Joy Axivane Fans, mounted in series aft of the 
test section. Each has an internal 1750 rpm, 100 horsepower motor. The fan blades are 
adjustable from 25 to 55 degrees. A set of multileaf-type inlet vanes immediately in front 
of each fan assembly allows the operator to control the mass flow rate and pre-swirls the 
flow in the direction of fan rotation, thereby reducing the power required to operate the 
wind tunnel. Flows of 10 to 250 feet per second can be produced. Figure 2 is an overall 
view of the wind tunnel. 
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2. Test Section 


The test section is 223 inches long and has a twenty-four inch square cross 
section. The top and bottom panels are constructed of continuous two-inch thick 
aluminum slabs. Access pons through the bottom slab permit the mounting of hotwires 
and pitot-static probes; they may be sealed off if not in use. The front sidewall is 
constructed of three two-inch thick stress-relieved Lucite panels, each secured by twelve 
bolts to the top and bottom slabs and to the vertical frames that divide the test section 
into three sections. The three front wall panels can be raised by hydraulic actuators to 
allow access to the test section. The rear side panels are two-inch thick plywood, and can 
be used for mounting internal or external instrumentation. The extraordinarily heavy 
construction of the test section is intended to minimize deflection due to rapidly changing 
static pressure in an oscillating flow. 

3. Rotating Shutters 

The method used to generate the oscillating flow was first used by Karlsson 
[Ref. 1], and later uy Miller [Ref. 2], The scheme consists of four equally spaced 
shutter vanes, mounted on horizontal shafts which span the wind tunnel horizontally just 
aft of the test section and protrude through the rear vertical side. The bottom shaft is 
driven by an electric motor, while the other shafts are each driven by the one immediately 
adjacent by means of a timing belt and pulley arrangement. This serves to keep all four 
vanes in phase as they rotate. Alternatively, the vanes can be locked in the horizontal 
plane to allow the tunnel to produce a s'eady flow. The electric motor’s speed can be 
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varied by the operator to produce an oscillation at any frequency from 2 to 240 Hz. The 
amplitude of flow oscillation about the mean can be altered by varying the width of the 
shutter vanes, and thus the degree of test section closure. In this work, shutter vanes of 
six-inch width were used, resulting in nearly 100 percent closure. With a mean flow 
velocity of 50 feet per second, oscillation amplitudes of approximately 35 feet per second 
RMS about the mean velocity were produced. Figure 3 is a photograph of the shutter 
vanes. 

B. MODEL 

The airfoil model was identical to that employed by Banning [Ref. 3]. It 
is a NACA 63-010 section modified by a straight line fairing from 60 percent chord to 
th trailing edge. The model has a 24-inch span and a 6-inch chord. Twenty-three 
pressure taps are distributed chordwise at mid-span on the top surface of the airfoil. Two 
additional taps are located on the underside, directly opposite ports number 11 and 16, 
to facilitate the calibration of a 0° angle of attack. Figure 4 shows the general shape of 
the airfoil and the location of the pressure ports. 

The airfoil model was pivoted at mid-chord for support and angle of attack control. 
Angle of attack was adjusted by means of a vernier mechanism, permitting the operator 
to set any desired angle from -40 to +40 degrees relative to the tunnel centerline. The 
pressure port leads penetrated the rear side wall through the angle of attack mechanism. 
They were then connected to the Scanivalve pressure measurement equipment, which is 
described below. Figure 5 is a photograph of the model mounted in the test section. 
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C. INSTRUMENTATION 


1. Velocity Measurement 


Two methods were used to measure flow velocity during experimental runs. 
First, a pitot-static tube was installed approximately 15 inches ahead of and 6 inches 
below the leading edge of the model, 4 inches from the side wall. 1 *.e pitot-static tube 
was connected to a water micro-manometer and was used to establish a known steady 












flow speed in order to calibrate the hotwire anemometer, and for measuring the mean 
flow velocity in oscillating flows. It should be stated that even with oscillations at 10 Hz, 
which produced velocities ranging from 0 to 80 feet per second, the mean velocity of 50 
feet per second was displayed steadily by the water manometer. All pressure fluctuations 
due to the velocity oscillation were damped out by the pitot-static tube and manometer. 

A hotwire anemometer was employed to measure the instantaneous velocity 
during all experimental runs. The hotwire was installed approximately 3 inches ahead of 
and 8 inches below the model leading edge, 6 inches from the side wall. The hotwire and 
pitot-static tube are visible in Figure 5. 

The linearized, constant temperature hotwire was described by Miller 

[Ref. 4], 

2. Pressure Measurement 

Instantaneous pressure was measured by means of a Scanivalve pressure 
transducer and signal conditioner. The Scanivalve unit was mounted on a free-standing 
mobile dolly and positioned adjacent to the model’s pressure port tubes where they exited 
the test section. This arrangement served to isolate the Scanivalve from the wind tunnel, 
thereby minimizing mechanical coupling and signal noise due to wind tunnel vibration. 

The model’s twenty-five pressure ports were connected to the Scanivalve by 
a set of uniform length thin-walled stainless steel tubes, each 24 inches long and with 
0.047 inches inner diameter. An additional lead was connected to a spare Scanivalve port 
for use in calibrating the pressure transducer. 
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The uniform length of the model’s pressure port tubing allows a single transfer 
function to be applied to the instantaneous pressure measurements from any port to 
account for phase lag and amplitude damping in the pressure signal due to the pressure 
port plumbing. A calibration was performed on the tubing to determine the transfer 
function necessary to convert the pressure signals received at the Scanivalve back to a 
signal representative of the pressure sensed at the model. An acoustic driver was used 
to generate a pressure signal of varying frequency at one end of the tubing. The pressure 
was measured at the source and at the remote end of the tube, and the resulting amplitude 
ratio and phase shift were plotted as a function of frequency. This technique has been 
described in detail by Banning [Ref. 3]. The results of this calibration are depicted in 
Figures 6 and 7. 

3. Data Acquisition System 

Instantaneous pressure and velocity measurements were recorded and stored 
on separate 64-K capacity recording boards installed in a 12 MHz 80286-based 
microcomputer. The data acquisition program was cued by an electrical trigger, which 
was actuated mechanically by the rotating shutter vanes. The data acquisition software 
armed the two data recording boards and, when triggered, the boards simultaneously 
recorded 1024 samples of pressure and velocity over a period of 1.024 seconds. This was 
programmed to match the oscillation frequency of 10 Hz, which was selected by the 
operator at the beginning of the experimental session. 

The pressure and velocity signals were output continuously from the 
Scanivalve pressure transducer and the hotwire anemometer, respectively, as DC voltages 
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corresponding to a percent deflection of full scale, where full scale refers to a full scale 
range on the data acquisition boards. In both cases, full scale deflection was limited to 
plus or minus five volts. As mentioned, the pressure and velocity analog signals were 
sampled at slightly greater than 100 times the 10 Hz frequency of the oscillating flow to 
produce a sample size of 2 N , which was necessary for compatability of the data 
acquisition system. The recording boards are capable of data acquisition up to a data rate 
of 1 MHz. 

The recording boards are of 12-bit resolution, capable of differentiating 1 part 
in 4096. By calibrating the pressure signal recording board response for a maximum 
signal of plus or minus 1 psi, an ability to resolve 0.0005 psi was achieved. The velocity 
recording board was calibrated to a maximum velocity of plus or minus 125 feet per 
second, therefore its resolution was 0.061 feet per second. 

4. Controls 

The operator’s console consisted of several panels of instrumentation which 
enabled the experimenter to control and monitor nearly every aspect of the experiment 
in progress. Figure 8 is a photograph of the console system. On the far left is a 
Tektronix Type 555 Dual Beam Oscilloscope, used to monitor the pressure and velocity 
signals in real-time. On the left end of the left console is the control panel for a the ten- 
channel linearized hotwire system. Only one channel was used in the present work. To 
the right of the hotwire system is a Singer Model 303-01 RMS Voltmeter, used to 
monitor the RMS velocity fluctuation during oscillating flow experiments. The DC 
voltmeter to its right may be used to monitor the mean velocity. On the upper right side 
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of the left console is a Dynasciences Signal Counter, which continually measured and 
displayed the rotating shutter frequency. Below that is an Analog Devices AD2501 
digital thermometer, which displayed freestream temperature in degrees Fahrenheit, 
measured by a thermocouple just downstream of the model. On the lower panel are 
located the hydraulic switches used to raise and lower the test section access doors, the 
mass flow rate adjustment sw itches which control the degree of closure of the multileaf 
vanes, the shutter activation and speed controls, i the emergency shutoff switches for 
the two electric fan motors. Tht ater manometer used to establish tunnel velocity for 
calibration and to measure freestream mean velocity in oscillating flows sits atop the desk 
in front of the left console. 

The right-hand console incorporates the power supply for the data acquisition 
trigger, the Scanivalve Position Display and solenoid controller, and the data 
acquisition/reduction microcomputer. The computer is an Everex System 1700 Model-C 
AT clone with an Intel 80286-12 CPU chip and an 80287 math coprocessor chip. Also 
seen are an oscilloscope and an Ampex PR2200 magnetic tape recorder, w'hich may be 
used to monitor and record analog signals. On the table next to the console tire an HP 
7470A plotter, IBM Proprinter II and an IBM PC/XT microcomputer, used for additional 
data reduction capability. 
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m. EXPERIMENTAL PROCEDURES 


A. CALIBRATION 

1. Flow Velocity 

Prior to each run, the hotwire anemometer was calibrated sing a water 
manometer to establish a known velocity for comparison. Calibration points were taken 
at 0 and 100 feet per second for every run. First, a zero velocity was assured by placing 
a cover over the hotwire probe. In this condition, the hotwire current was adjusted at 
maximum gain to obtain an output of 0 volts, to an accuracy of 0.0001 volts. Then the 
wind tunnel was started and the flow adjusted to generate a water column of the proper 
height for 100 feet per second, based on local barometric pressure and temperature. The 
hotwire gain control was adjusted to obtain an output of 4.0000 volts to an accuracy of 
0.0005 volts. To avoid over-voltage of the data acquisition system, a limit of 5 volts, 
corresponding to 125 feet per second, was established for velocity measurements. For 
calibration, these voltages were determined by the digital voltmeter in the linearized 
hotwire circuit. 

2. Pressure Measurement 

Pressure measurement calibration was performed by first setting the Scanivalve 
pressure transducer output equal to zero with the calibration pressure port open to the 
atmosphere. Then, the pressure transducer voltage gain w'as adjusted to read -2.50 volts 
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with an external pressure input of 0.50 psi above atmospheric pressure applied, as read 
on a water column manometer. Thus, the voltage output of the Scanivalve would 
represent: 

p =p~p (It 

Scanivalve * Model Atmosphere ' ' 

The maximum voltage span for the Scanivalve was set at plus or minus 5 volts, allowing 
a pressure equal to atmospheric plus or minus IT* psi to be measured without causing 
over-voltage of the data acquisition system The pressure calibration voltages were 
determined by using a Fluke Digital Voltmeter, with an accuracy of 0.001 volts. 

3. Data Acquisition System 

The data acquisition system boards were calibrated for a true zero reading at 
a 0.0 volt input, and set to accept a 10-volt span, corresponding to a reading of 2048 
digital counts at +5.00 volts ana -2047 digital counts at -5.00 volts, for both pressure and 
velocity measurements. The calibration results are presented in Figure 9, which shows 
a linear, symmetrical response centered within 0.1 percent of true zero. 

B. TUNNEL OPERATION 
1. General 

Once calibration procedures were completed, the tunnel velocity, oscillation 
frequency and model angle of attack were set at the desired values. The operator supplied 
information regarding tunnel operating conditions and ambient conditions (pressure, 
temperature, etc.) to the data acquisition program for inclusion in the data file header. 
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Figure 9. Data Acquisition Boards Calibration Results. 

The computer controlled the actual data acquisition process, pausing after reading the first 
23 pressure pons to allow the operator to rotate the model to a negative angle of attack. 
Thus, the same pressure taps represented both the upper an lower surface of the 
symmetrical airfoil. The tunnel velocity and Scanivalve pressure and port selection were 
monitored continually by the operator to ensure proper operation throughout the run. 

2. Steady Flow 

For steady flow operation, the rotating vanes were secred fully open in the 
houzontal position. A flow velocity of 100 feet per second was established as measured 
by the pitot-static tube and water manometer, with the model at the desired angle of 
attack. Model angles of attack from 0 to 20 degrees, in increments of 5 degrees, were 
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investigated to establish baseline values for aerodynamic coefficients and steady-state 
pressure distribution as a function of angle of attack. The chord Reynolds number for 
steady flow work was 306.000, where the Reynolds number was given by: 

R = ( 2 ) 

M- 

3. Oscillating Flow 

All oscillating flow runs were made at an oscillation frequency of 10 Hz. The 
mean flow velocity was established at 50 feet per second as measured by the pitot-static 
tube and water manometer, with the vanes rotating and the model fixed at the desired 
angle of attack. The velocity oscillations were of the order of plus or minus 70 percent 
of the mean value. Model angles of attack from 0 to 35 degrees, in increments of 5 
degrees, were investigated. A Reynolds number of 152,000 and a Strouhal number of 
0.314, referenced to (he mean flow velocity and the model chord length, were used for 
all oscillating flows. Strouhal number was determined by the relation: 


S - 



U 


(3) 


A typical Strouhal number for rotary wing aircraft is on the order or 0.025. 
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IV. DATA ANALYSIS 


A. DATA ACQUISITION 

An automated data acquisition program controlled the progress of the experimental 
runs under the supervision of the operator. Because the experimental data, consisting of 
the time-dependent pressures and velocities, were stored in order of occurrence as 
voltages on the data acquisition boards, a data acquisition program was written in Quick 
BASIC, which permitted the interrogation of specific memory locations on the recording 
boards and retrieval of discrete pieces of information related to the control of the 
experiment. 

During the data acquisition cycle, data was taken for 1.024 seconds at a time, 
which spans 10 cycles of flow oscillation. The signals were ensemble averaged to 
attenuate the effect of transients in the flow, and the time-dependent portion of each 
signal was separated from its respective mean value for each pressure port location. The 
time-dependent portion of the pressure signal then underwent a Fourier transform and was 
divided by the previously determined transfer function, correcting for the amplitude 
damping and phase shift associated with the pressure tap connecting tubes. Next, an 
inverse Fourier transform was applied. This resulted in a discrete representation of the 
pressure distribution at each port and the corresponding velocity profile over a single 
cycle of flow oscillation. By separating the time-dependent portion of the two signa. s, 
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greater accuracy could be achieved in the Fourier transforms and transfer function 
processes. 

Output files were formatted such that the header information was listed at the 
top of the file, and could be accessed for data reduction. This was followed by the 46 
average pressures, 46 average velocities, 4600 values of the time-dependent pressure 
signal and 4600 values of the time-dependent velocity signal. A naming convention was 
adopted for easy file identification, wherein the output files were named RUNxyy.DAT, 
where ’x’ refers to the run series (0 for steady flow, 1 for oscillating flow) and ’yy’ 
refers to the angle of attack, in degrees. Figure 10 is a block diagram of the data 
acquisition program, and a complete program listing appears in Appendix A. 

B. DATA REDUCTION 

The data reduction program was written in FORTRAN for simplicity and for speed 
in repetitive calculations. The program used the RUNxyy.DAT file in conjunction with 
another file, XOVRC.DAT, containing the pressure port chordwise locations, and wrote 
information to nine output files including the upper and lower surface pressure, velocity 
and pressure coefficient as a function of chordwise position and cycle time, and regarding 
the C L , C D and C M c/4 values as a function of cycle time and angle of attack. The pressure 
and velocity files were examined to ensure correlation with what was observed on the 
oscilloscope during the experiments. The C P , C L , C D and C M c« files are included in this 
report as two-dimensional curves for steady flow results and as three-dimensional surface 
plots for oscillating flow results, in Chapter V. 


22 






READ TRANSFER FUNCTION 

This correlates the pressure 
read by the Scanivaivc to the 
actual pressure, P(t), which is 
experienced at the model 


READ HEADER INFORMATION 

Program queries operator for 
run number ’xyy\ AOA, P ABI1 , T 
and flow oscillation frequency 


SET SCANIVALVE TO HOME 


ADVANCE SCANIVALVE TO PORT 1 


INITIALIZE RECORDING BOARDS 

AND ANALOG-TO-DIGITAL SYSTEM 

Performed once, prior to each 
data run, to establish various 
recording board control parameters 

FOR FORT = ! TO 23 

READ PORT NUMBER 


RECORD PRESSURE & VELOCITY 

1024 values of P and V are 
recorded simultaneously by 
die recording boards. 

NEXT FORT 

LOAD P, V INTO ARRAYS AND PROCESS 

Ensemble average 10 cycles’ worth 
of P, V values, separate P s . iAN and 

Vmf-an from P(t) and V(t). Apply FFT, 
transfer function & IFFT to pressure 
signal only. Scale P 'o !bf/'ft : , V to 
ft/sec. 


PAUSE TO ALLOW OPERATOR TO SET 

NEGATIVE AOA, REPEAT LOOP TO 

RECORD LOWER SIDE DATA 


WRITE OUTPUT FILE RUNxyy.DAT 


Figure 10. Data Acquisition Program Flow Chart. 
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As previously described, the Scanivalve output was set to 


Pscanl valvG~ PMode 1 Pa tm 


( 1 ) 


Total pressure was found to be equal to atmospheric pressure, due to the low inlet 
velocity through the high solidity inlet screens. Thus, 

P 0 =PAUn=P-+Q- < 4 > 


The coefficient of pressure is given by 


r, _ PModel Pm 
L- n 


(5) 


Substituting (1) and (4) into (5) yields 




t-^scan + -^ , .Atm 


1 - [P Acm -gJ 




( 6 ) 


which further simplifies to 


C p - 


_ PSeem + _ Pscan 


Qm 


+ 1 


(7) 


This result was applied to each pressure value to compute the corresponding C P value. 
The resulting distribution with respect to cycle time and chordwise parameter was 
numerically integrated in the following manner 
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c« 



( 8 ) 


to obudn a normal forcv coefficient, C N . This was resolved into lift and drag components 
by applying the cosine or sine of the angle of attack, as appropriate. Likewise, the 
following computation 



f lC p -C p ] * [ 0 . 25 - — ] d— 
J 1 « C C 



c 


( 9 ) 



was applied to obtain the moment coefficient about the quarter chord, Cmc/ 4. A block 
diagram of the data reduction program appears in Figure 11, and a complete program 
listing appears in Appendix B. 

C. UNCERTAINTY ANALYSIS 

As previously stated, the data acquisition system was calibrated to determine the 
linearity of response to a known input and the accuracy of the digital count with reference 
to true zero. The calibration data is reported in Figure (9). True zero appeared at -5 of 
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2047 counts on the velocity board, and at -3 of 2047 counts on the pressure board. The 
slope of both boards was linear and consistent with the ideal slope, corresponding to a 
count of -2047 at an input of -5.0 volts and a count of +2048 at an input of +5.0 volts, 
to within 0.1 percent. With the range of the pressure board set at plus or minus 1.0 psi, 
or 144.0 psf, this corresponds to a maximum error of 0.144 psf, or 0.001 psi in reading 
a full-scale deflection pressure of 1.0 psi. Similarly, for a full scale velocity reading of 
125 feet per second, the maximum error would be 0.125 feet per second. Since the 
maximum values encountered were approximately 0.3 psi and 85 feet per second, the 
uncertainty in these measurements is of 
the order of 0.1 percent. 

The sensitivity of the Scanivalve 
pressure transducer is guaranteed to be 
17.5 millivolts, plus or minus 3 millivolts, 
as published in the Scanivalve calibration 
certificate. Given the voltage and 
pressure scaie in use, this results in a 
contribution to the uncertainty of 0.41 
percent due to the Scanivalve. 

The water manometer is marked in 
graduations that permit the operator to 
distinguish 0.005 centimeters of water. 


OPEN RUNxyy.DAT, XOVRC.DAT 
READ NECESSARY HEADER INFO 


RECONSTRUCT PRESSURE, VELOCITY 
DISTRIBUTION FROM P MEAN , V^, 
P(t) & V(t) 


COMPUTE CpUPPER AND C P LOWER 
AND STORE IN ARRAYS 


WRITE P, V, C P OUTPUT FILES 
FOR THIS RUN NUMBER 


INTEGRATE C P DISTRIBUTION 
TO OBTAIN C N , CmC/4 


COMPUTE C L , C D AND WRITE 
THESE PLUS CyC/4 TO THEIR 
RESPECTIVE OUTPUT FILES 


Figure 10. Data Reduction Program Flow 
Chart. 
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For an oscillating flow of 50 feet per second mean velocity, this corresponds to an 
uncertainty of 0.41 percent. 

The cumulative result of these uncertainties in calculation of C P , assuming the worst 
case of an overestimation of the Scanivalve pressure reading combined with an 
underestimation of the dynamic pressure, may be calculated starting with the Scanivalve 
uncertainty: 

^Scanivaive = l • ®1 *Pxiue (H) 

Then the effect of recording board nonlinearity may be written: 

■Pj?ec.Board~l • • 0051*P rrue (12) 

The effect of uncertainties in the value of the dynamic pressure, q, may be written: 

Qrnpui'® • 9^66 *Q Trug ( 13 ) 

These results may then be used to calculate the uncertainty in C P : 

C -C 

Uncertainty= - w/errors —(14) 

p Tiue 


or: 


• Roc.Board 
QInput 


1] ~ [ 


■ True 


+ 1 ] 


i Tiue 


1] 

^Tzue 


(15) 
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Finally, using equations (12), (13) and (15), one may write: 


Uncertainty = 


[ 


1 ■ 00 51*P True 
0 .9966*q True 


+!]-[ 


Qtzuo 


1] 

Qti uo 


( 16 ) 


or simplifying: 


[0.0085* F ^ ue ] 

Uncer tain ty- - ?— ue . (17) 

[J"H£ +1 ] 

Qtiub 


Thus the uncertainty in the calculation of C P approaches a maximum value of 0.85 percent 

P 

as the ratio — >1 . 

Q 


The uncertainty in the coefficients, C L , C D and C M c/4, which are derived by 
integrating C P , may be expected to be of a similar value. All computation in the data 
reduction program was performed using FORTRAN double precision, to avoid 
propagation of roundoff errors. 
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V. RESULTS AND DISCUSSION 


A. STEADY FLOW RESULTS 

The C P distributions for the five steady flow runs are presented below in Figures 
12-16. The plots in these figures were generated by averaging the corresponding C P 
output files with respect to cycle time, t/T, since there were no oscillations in the flow. 
The resulting data points are all displayed, connected by a cubic spline curve of medium 
tension. 



Figure 12. Pressure Coefficient Distribution at 0° Angle of Attack in Steady Flow. 
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The pressure distributions in steady flow are as one might expect for a symmetrical 
airfoil of this type. It can be seen that a slight leading edge separation bubble exists 
between ports 2 and 4, indicated by irregularities in several of the C P curves at that 
location. The mean pressure distributions at 5, 10 and 15 degrees angle of attack are 
quite similar to those previously reported by Banning [Ref. 3|. C P distributions at 0 and 
20 degrees angle of attack were not reported by Banning. 

The steady flow performance coefficients, C L , C D and C M c/4, are shown below as a 
function of angle of attack in Figures 17-19. The lift coefficient curve agrees well with 
that reported by Banning, who did not report drag or moment coefficients. The lift 
coefficient curve exhibits classical behavior for steady flow, with stall occurring at 



Figure 17. Lift Coefficient in Steady Flow. 
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approximately 12 degrees. The offset from the origin at 0° angle of attack is attributed 
to possible pressure tap errors near the leading edge. 


B. OSCILLATING FLOW RESULTS 

The C P distributions a function of dimensionless time, t/T, and chordwise position, 
x/c, for angles of attack from 0 to 35 degrees, and their respective velocity profiles, are 
presented in Figures 20-43. The C P distribution surface plots were generated by first 
creating a grid consisting of chordwise data at intervals of x/c = 0.025 and timewise data 
at intervals of t/T = 0.01. Thus some interpolation had to be made to create a 3900-point 
grid out of the 2300 existing data points. This was done using Surfer, a three- 
dimensional surface plot software package, which uses an inverse distance squared routine 
wherein 


Z= 



( 16 ) 


In this equation, Z is the C P value of the grid point to be calculated, d is the distance to 
the a neighboring grid point, and n is the number of nearest neighbors to be used in the 
calculations. A grid size of 0.025 by 0.02 was used, to ensure that only the existing data 
points were used, to the maximum extent possible, with a minimum number of grid points 
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being interpolated to fill in between any two actual data points. The surface plots are 
presented with only chordwise lines, that is, lines of constant time t/T. Additionally, the 
plots have been rotated and elevated so as to give the best view of the major features. 
The reader should therefore take care to note the orientation of the axis scales. 

As is readily apparent, the pressure distribution over the airfoil was primarily a 
result of the oscillating flow and its concurrent oscillation in static pressure. As a 
consequence, chordwise pressure variations are very small compared to these induced 
static pressures. However, in spite of the lack of an orderly chordwise variation in C P , 
one can see by comparing the maximum and minimum values of -C P found on the upper 
and lower surfaces that the upper surface appears to have higher peaks -C P s at t/T = 0.13 
and 0.55, whereas deeper valleys are found at t/T = 0.40 and 0.80 in the lower surface 
C P surface plots. The net result, when the curves are integrated, is a non-zero normal 
force coefficient. 

Some of the distributions exhibit random anomalies at a single chordwise location, 
usually just aft of mid-chord. These anomalies are most noticeable when they exceed the 
peak value of negative C P and when they bridge the early, more shallow valley. In terms 
of the velocity cycle, this corresponds to the region in the velocity profile at which the 
flow deceleration suddenly becomes quite steep and quickly reaches a minimum. This 
feature of the C P surface plots is attributed to the onset of separation. The condition 
remains generally only until the flow is re-established in a steep acceleration. This 
phenomenon is most noticeable on the upper surface of the airfoil, as would be expected. 


Figure 20. Upper Surface Pressure Coefficient at 0° Angle of Attack in Oscillating 
Flow. 



Figure 21. Lower Surface Pressure Coefficient at 0° Angle of Attack in Oscillating 
Flow. 
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Figure 23. Upper Surface Pressure coefficient at 5° Angle Of Attack in Oscillating 
Flow. 



Figure 24. Lower Surface Pressure Coefficient at 5° Angle of Attack in Oscillating 
Flow. 
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Figure 28. Velocity Profile, RUN110. 
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Figure 29. Upper Surface Pressure Coefficient at 15° Angle of Attack in Oscillating 
How. 



Figure 30. Lower Surface Pressure Coefficient at 15° Angle of Attack in Oscillating 
Flow. 
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Figure 32. Upper Surface Presure Coefficient at 20° Angle of Attack in Oscillating 
Flow. 



Figure 33. Lower Surface Pressure Coefficient at 20° Angle of Attack in Oscillating 
Flow. 
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Figure 35. Upper Surface Pressure Coefficient at 25° Angle of Attack in Oscillating 
Flow. 



Figure 36. Lower Surface Pressure Coefficient at 25° Anglp of Attack in Oscillating 
Flow. 
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Figure 38. Upper Surface Pressure Coefficient at 30° Angle of Attack in Oscillating 
Flow. 



Figure 39. Lower Surface Pressure Coefficient at 30° Angle of Attack in Oscillating 
Flow. 
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Figure 41. Upper Surface Pressure Coefficient at 35° Angle of Attack in Oscillating 
Flow. 



Figure 42. Lower Surface Pressure Coefficient at 35° Angle of Attack in Oscillating 
Flow. 
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It is most common in the lower angles of attack, particularly from 5 to 15 degrees, and 
less so in the higher angles of attack, until its reappearance at 35 degrees. This is 
consistent with the known phenomenon of increased C L at high angles of attack, which 
is thought to be caused by delayed onset of dynamic stall in oscillating flow. 

The oscillating flow performance coefficients are plotted as a function of angle of 
attack and cycle time in Figures 44 - 49. These plots were also generated with Surfer, 
but with an initial grid size of 8 by 25, which used every fourth data point in the 
timewise direction. This grid was then spline smoothed, inserting two grid points in 
between each existing data point in both the angle-of-attack-wise and timewise directions. 
The lines in the first three surface plots are constant with respect to the cycle time, and 
the lines in the last three are constant with respect to angle of attack. The most obvious 
feature of all of the performance coefficient curves is their apparently random, 
unpredictable appearance. In spite of the closeness with which the coefficient of pressure 
adhered to the dominant freestream static pressure oscillations, there is no observable 
pattern that repeats in time at any constant angle of attack. Similarly, there is no 
evidence of steady or quasi-steady flow behavior. This observation can be verified by 
following any constant-time line and noting that it does not trace out a reasonable 
facsimile of a steady flow lift, drag or moment coefficient curve. This result is not 
surprising in light of the absence of any traditional chordwise pressure variation. 

To a small degree, there is a timewise variation in the lift coefficient similar to that 
in the pressure coefficient plots. In general, the value of the coefficient rises and falls 
with respect to time in a manner similar to C P , forming two ridges and two valleys. This 
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Figure 44. Lift Coefficient in Oscillating Flow (Constant 
Time Lines) . 
















Figure 46. Quarter Chord Moment Coefficient in Oscillating 
Flow (Constant Time Lines) . 
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Figure 49. Quarter Chord Moment Coefficient in Oscillating 
Flow (Constant AOA Lines). 
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pattern is somewhat distorted throughout the region between 5 and 15 degrees angle of 
attack, which is consistent with the separations noted previously in the pressure coefficient 
plots. Additionally, there is a general trend of increasing lift coefficient from 0 to 10 
degrees, followed by a lower, flat area at 15 to 20 degrees, and a subsequent increasing 
trend out to 35 degrees angle of attack. 

It is apparent that the drag coefficient increases in an exponential fashion with 
increased angle of attack. Likewise, the moment about the quarter chord appears to 
become substantially more negative at higher angles of attack. The drag and moment 
coefficient curves are more difficult to interpret in general. 

The mean values of the lift, drag and moment coefficients are presented as a 
function of angle of attack only in Figures 50 - 52. These graphs are also plotted such 



Figure 50. Mean Lift Coefficient in 10 Hz Oscillating Flow. 
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Figure 51. Mean Drag Coefficient in Oscillating Flow. 



Figure 52. Mean Quarter Chord Moment Coefficient in 
Oscillating Flow. 
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that the original data points are displayed, each representing the arithmetic mean value 
of the respective coefficient at constant angle of attack, connected by a cubic spline curve. 
These performance coefficients exhibit the well known ’’average” behavior of a 
symmetrical airfoil in oscillating flow. In particular, the lift coefficient appears to 
increase at angles of attack in excess of the steady flow stall angle, even out to 35 
degrees. The magnitude of the lift coefficient at lower angles of attack is slightly lower 
than the corresponding values in steady flow, however the peak value of 2.9 at 35 degrees 
is well in excess of the steady flow peak C L of 0.78, found at 10 degrees angle of attack. 
The oscillating flow mean drag coefficient is increased at every angle of attack over the 
corresponding steady flow drag coefficient, however the general trend of the two plots is 
to increase with increasing angle of attack. Likewise, the oscillating flow mean moment 
coefficient curve is similar in shape to that of the steady flow moment coefficient but 
with exaggerated values at the higher angles of attack. 
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VI. CONCLUSIONS 


The steady flow pressure and performance coefficient results lend credence to the 
experimental procedures, data acquisition and data reduction in this work. The pressure 
coefficient curves show entirely predictable chordwise variation, which is well known to 
produce lift on an airfoil. Conversely, the oscillating flow pressure coefficients exhibit 
no such traditional chordwise variation, and in fact indicate that the instantaneous pressure 
distribution on a symmetrical airfoil in an oscillating flow is dominated by the freestream 
static pressure. The noted exceptions to this observation are the occasional flow 
separations at the lower angles of attack, analagous to the onset of separation in steady 
flow. 

The oscillating flow results also indicate that the aerodynamic performance of a 
symmetrical airfoil in an oscillating flow is in no way a quasi-steady phenomenon. This 
is evidenced by the fact that the surface plots of the lift, drag and moment coefficients 
have a very random appearance, in spite of the closeness with which the pressure 
coefficient follows the freestream static pressure. It is not intuitively obvious, based on 
the shape of the velocity profile and the pressure coefficient surface plots, that the 
performance coefficients would behave as they do. 

This randomness notwithstanding, the mean values of the same performance 
coefficients attest to the known behavior of increased lift coefficients in oscillating flow. 
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As previously noted, the general trends of the lift, drag and moment coefficients in 
oscillating flow resemble their steady flow counterparts, albeit with exaggerated values. 
One must conclude, however, that the instantaneous values of these coefficients cannot 
be assumed to be a continuum of quasi-steady states. 
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APPENDIX A 


The following is the data acquisition program, TKDATA.BAS, which is written in 
Quick BASIC. 


’PROGRAM TKDATA This program contains SUBROUTINES to control the SCAN1VALVE, 
’setup the SRA 1200 RECORDER boards, take dat, download that 
’data into arrays VRAW, representing the raw velocity data, 

’PRAW, representing the raw pressure data, process that data, 

’and write the output to file and/or plot it on the screen. 

’These SUBROUTINES and their respective functions are: 

’HOME - Returns SCANIVALVE to HOME 

’HEADER - Inputs data for a file header 

’PORTNUM - Returns SCANIVALVE number (PORTNO) 

TNIT - Initializes recorder boards 

’RECORD - Records data 

’getstanadd - Returns starting address of record 

’getdai - Returns data samples 

’PLOT - Plots data on screen 

’LOAD - Loads two arrays, VRAW and PRAW, (in pseudo- 
’complex format) 

’MEAN - Computes the mean of the arrays PRAW and VRAW, 

’subtracts the means and scales the values 
’FFTPRAW - Replaces PRAW with the FFT of PRAW 
’IFFTPRAW - Replaces PRAW with the IFF! of PRAW 
’PORT - Advances the SCANIVALVE port 
'RiRTHEi A - Converts an array from real-imaginary pairs 
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’to r-theta pairs 

’RTHETAR1 - Converts an array from r-theta pairs to 
’real-imaginary pairs 

’DIVIDE - Replaces PRAW with PRAW/XFER 
’XFERINPUT - Inputs the transfer function array XFER 
’from an external file XFER.FCN 
’NEGALPHA - Halts data collection to allow resetting 


’angle of attack to MINUS ALPHA 
’ARRAY - Loads mean values of pressure and velocity 
’into vectors PMEAN and VMEAN and pressure 
’and velocity data into arrays P and V 
’DISK - Outputs pressure, velocity and header data 

’to a file called DATA.XXX where XXX is the 
’Run Number. 

VSCALE =125! 

PSCALE = 144! 

OPTION BASE 1 ’Sets lower bound of subscripts = 1 

DIM physicaladdrcss(2) ’Physical adresses of the recorders 

DIM inputscale(2) ’Full scale input voltage for each 

’channel. May be 1, 2, 5, or 10 volts 


DIM starts(2) 

DIM VRAW(2048) 
DIM PRAW(2048) 
DIM XFER(2048) 
DIM PMEAN(46) 

DIM VMEAN(46) 
DIM P(46, KXfi 
DIM 100) 
GOSUB XFERINPUT 
GOSUB HEADER 
GOSUB HOME 
GOSUB PORT 
GOSUB INIT 


’Array containing starting addresses 

’Array containing raw velocity data 
’Array containing raw pressure data 
'Array containing the transfer fnen. 
'Vector containing mean pressure 
’Vector containing mean velocity 
’Array containing pressure data 
’Array containing velocity data 

’Inputs transfer function XFER 
’Calls for header information 
’Sets SCANIVALVE to PORT 0 
'Advances SCANIVALVE to PORT 1 
'Use once before using recorders 


7 IlllllllllllllllltlllllllllllllllllllllllllllllllilllllllllllllllllllllllH 


’RECORDER CONTROL PARAMETERS — Set for specific application - 
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These parameters are: 

rate, the sample rate in mico-seconds, from 1 to 32767 
samples, the number of samples per record, from 1 to memmax 
psamples, die number of samples that preceed the trigger from 
1 to memmax, but less than samples. 
modeS, trigger slope control, either "PLUS" or “MINUS" 
setpoint, for trigger threshold, range is -100 to 100 % 

optional: 

inputscale(n) controls programmable gain for channel n. 

Valid with values of 1, 2, 5 and 10; n may be either 1 or 2. 


SAMPLES = 1024 

’Number of samples per record, (must 


’be of the form 2**n for processing) 

psamples = 1000 

’Number of pretrigger samples 

setpoint = 10 

’Trigger setpoint in % of full scale 

modeS = "PLUS" 

’Trigger modeS, (either "PLUS" 


’or "MINUS") 

rate = 1000 

’Sample rate in u-seconds per point 

inpuiscale(l) = 5 

’full scale for channel 1 is 5 volts 

inputscale(2) = 5 

'full scale for channel 2 is 5 volts 


’ llllllllllllllllllllllllllllllllllllilllllllllllHIIIIIIIIIIIIIIIIIIIIIIIIil 


Begin first data acquisiti i 

loop through top 23 ports: 

FOR J = 1 TO 20000 

'Delays data acquisition to allow 

NEXT J 

’SCANiVALVE to settle 

FOR P = 1 TO 2. 


GOSUB PORTNUM 

’Reads current port number 

GOSUB RECORD 

’Records data for one port 

’GOSUB PLOT 

’Plots recorded data for debugging 

’END 

’Test point for debugging 

GOSUB LOAD 

’Loads data into arrays VRAW & PRAW 


’of length 2048 in pseudo complex 
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'form 


’END ’Test point for debugging 


GOSUB PORT ’Advances SCANT VALVE to next port 

GOSUB MEAN ’Computes the mean of each array, 

’subtracts the means (PBAR and VBAR), 

’and scales all values with PSCALE 
’and VSCALE so that output is in psf 
’and ft/sec. Arrays remain named PR AW 
’and VRAW 


’END 


’Test point for debugging 


GOSUB FFTPRAW 
GGSUo R1RTHETA 

GOSUB DIVIDE 
GOSUB RTHETARI 
GOSUB IFFTRAW 


’Replaces PRAW with its FFT 
’Transforms PRAW in pseudo-complex 

’form 

’Divides FFT of PRAW by XFER 
’Transforms PRAW back r-i pairs 
’Performs inverse FFT on PRAW 


’END 


’Test point for debugging 


GOSUB ARRAY ’Loads PBAR and VBAR into vectors 

’and PRAW and VRAW into arrays 


'END 


’Test point for debugging 


NEXT P 

’Next change angle of attack: 

GOSUB NEGALPHA ’Stops execution to change ALPHA 

’Then wc continue with data collection for the bottom 23 ports: 

GOSUB HOME ’Sends SCANIVALVE back to zero 

GOSUB PORT ’Advances SCANIVALVE to first port 

FOR J = ’ TO 20(KK) ’Delays data acquisition to allow 
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NEXT J 

FOR P = 1 TO 23 


’ SCANIVALVE to settle 


GOSUB PORTNUM ’Reads current port number 

PORTNO = PORTNO + 23 ’Adds 23 to PORTNO for no. 24 thru 46 

GOSUB RECORD ’Records data for one port 

GOSUB LOAD ’Loads data into arrays VRAW & PRAW 

’of length 2048 in pseudo complex 
'form 

GOSUB PORT ’Advances Scanivalve to next port 

GOSUB MEAN ’Computes the mean of each array, 

’subtracts the means (PBAR and VBAR), 

’and scales all values with PSCALE 
'and VSCALE so that output is in psf 
’and ft/sec. Arrays remain named PRAW 
’and VRAW 

GOSUB FFTPRAW ’Replaces PRAW with its FFT 

GOSUB RIRTHETA ’Transforms PRAW in pseudo-complex 

’form 

GOSUB DIVIDE Divides FFT of PRAW by XFER 

GOSUB RTHETARI ’Transforms PRAW back r-i pairs 

GOSUB IFFTRAW ’Performs inverse FFT on PRAW 

GOSUB ARRAY ’Loads PBAR and VBAR into vectors 

’and PRAW and VRAW into arrays 

NEXT P 

GOSUB DISK: ’Writes all header data, the vectors 

’PMEAN and VMEAN and the arrays P and 
’V to disk in a file RUNXXX.DAT 

END ’End of application code. 


’ SUBROUTINES WHICH MAY BE CALLED AS REQUIRED BY PROGRAM TKDATA 
INIT; ’Use this routine once at startup for variable and nardware init. 
k = 1024 ’k is 2 ** 10 

SAMPLES = 2000 ’init # of samples per record 

psamplcs = INT(SAMPLES / 4) ’init presamplcs at .25 total record 
setpoint = 50 ’trigger setpoint in % of full scale 
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modeS = "PLUS" 
rate = 4 

memmax = (64 * k) - 1 

tbc = &H205 

cwr = &H21E 

ctrO = &H206 

Ctrl = &H20E 

ctr2 = &H216 

clkie = 1 

FOR J = 1 TO 2 


’trigger modeS is either "PLUS" or "MINUS" or - 
’default sample rate in us. 

’recorder memory size 
’time base control register address 
’8253 control write register 
’8253 counter #1 
’8263 counter #2 
’8253 counter #3 

’for hardware hook not being used 


’set up list for physical address starting at &h200 
physicaladdress(J) = &H200 + (J * 32) - 32 
’Read board jumpers and set inputscale(j) 

’Absent boards get inputscale = 0 
’Others will return 1, 2, 5 or 10 

x = INP(physicaladdress(J) + 5) ’read jumpers w4, w5, w6, w7 
x = x AND &H70 


SELECT CASE x 
CASE 96 
inputscale(J) = 1 
CASE 80 
inputscale(J) = 2 
CASE 64 
inputscalc(J) = 5! 

CASE 48 
inputscale(J) = 10! 

CASE ELSE 
inputscale(J) = 0 
END SELECT 
NEXT J 

’Initialize time base and run through one cycle: 

rate = 1000 

’Toggle reset line 

OUT tbc, &HEF 

OUT tbc, &HE3 
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’Set counter 0 to mode 3 
OUT cwt, &H36 
’Set sample rate 
A = rate * 2 
ha = INT(A / 256) 
la = A - (ha * 256) 

OUT ctrO, la 

OUT ctrO, ha 

’Set counter 1 to mode 3 

’And divide by 4 

OUT cwr, &H76 

OUT Ctrl, 4 

OUT Ctrl, 0 

’Set counter 2 to mode 5 
OUT cwr, &HBA 
’Set number of samples 
B = (SAMPLES - psamples) / 4 
hb = INT(B / 256) 
lb = INT(B - (hb * 256)) 

OUT ctr2, lb 
OUT ctr2, hb 
’Set trigger level 

A = INT(setpoim * 1.28 * -1) + 127 -7 
OUT &H203, A 
’Arm time base in chosen mode 
’First create time base control word 
tbctw = &HF7 
’clkie always equals 1 

IF clkie = 1 THEN tbctw = tbctw AND &HFB 
’Test mode$( PLUS or MINUS ) 

IF modeS = "MINUS” THEN tbctw = tbctw AND &H7F 
OUT tbc, tbctw AND &HED 
’Wait loop to get pretrigger data 
ptdelay = psamples * rate * .000001 
A = TIMER + ptdelay + 1 
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’This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEYS = WEND 
OUT tbc, tbctw AND &HFE 
OUT tbc, tbctw 
’Force a trigger 

OUT tbc, tbctw OR 8 
FOR J = 1 TO 500: NEXT J 
OUT tbc, tbctw 

’This ends the time base first cycle and initialization.. 

RETURN 

RECORD: ’recording sequence control 
CLS 
’VIEW 

’WINDOW (0, 2048M5U, -2048) 

’LINE (50, 2000)-(461, -1850), lincol, B 
LOCATE 3, 27 

PRINT "RECORDING SEQUENCE CONTROL ” 

LOCATE 8, 25 

PRINT "1. ) Initialize Recorder Hardware " 

’First set the gain on boards 1 and 2. 

FOR J = 1 TO 2 
SELECT CASE inputscalc(J) 

CASE 1 
mask = &H9F 
CASE 2 
mask = &HDF 
CASE 5 
mask = &HBF 
CASE 10 
mask = &HFF 
CASE ELSE 
mask = &HFF 
END SELECT 
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IF J = 1 THEN maskl = mask 
OUT physicaladdress(J) + 5, mask 
NEXT J 

Toggle reset line 

OUT tbc, &HEF AND maskl 

OUT tbc, &HE3 AND maskl 

’Set counter 0 to mode 3 

OUT cwr, &H36 

’Set sample rate 

A = rate * 2 

ha = INT(A / 256) 

la = A - (ha * 256) 

OUT cuO, la 

CUT ctrt), ha 

’Set counter 1 to mode 3 

’And divide by 4 

OUT cwr, &H76 

OUT Ctrl, 4 

OUT Ctrl, 0 

’Set counter 2 to mode 5 
OUT cwr, &HBA 
’Set number of samples 
B = (SAMPLES - psamplcs) / 4 
hb = INT(B / 256) 
lb = INT(B - (hb * 256)) 

OUT ctr2, lb 
OUT ctr2, hb 
’Set trigger level 

A = INT(sctpoint * 1.28 * -1) + 127 - 7 
OUT &H203, A 
’Arm timebase in chosen mode 
’First create time base control word 
tbctw = &HF7 AND maskl 
’elkie always equals 1 

IF clkic = 1 THEN tbctw = tbctw AND &HFB 
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’Test mode$( PLUS or MINUS ) 

IF modeS = "MINUS" THEN tbctw = tbctw AND &H7F 

LOCATE 10, 25 

PRINT "2.) Pretrigger Delay " 

OUT tbc, tbctw AND &HED 
REM out tbc.tbctw 
’Wait loop to get pretrigger data 
ptdelay = psamples * rate * .000001 
A = TIMER + ptdelay + 1 

’This loop generates and error when timing through midnight 
WHILE (TIMER < A AND INKEYS = WEND 
LOCATE 12, 25 

PRINT "3.) Arm Recorder...wait for trigger" 

OUT tbc, tbctw AND &HFE 
OUT tbc, tbctw 
BS = "" 

A = 1 

WHILE (A o 0 AND INKEYS = "") 

A = INP(&H205) AND 3 

IF INKEYS o THEN 
OUT tbc, tbctw OR 8 
FOR I = 1 TO 500: NEXT J 
OUT tbc, tbctw 
END IF 
WEND 

LOCATE 14, 25 

PRINT "4. ) Triggered ...wait until recording complete" 

WHILE (A = 1 AND INKEYS = "") 

A = INPC&H205) AND 1 
WEND 

LOCATE 16, 25 

PRINT "6. ) Recording port number..."; PORTNO 
PRINT "6, ) Hit any Key to Exit" 

AS = "" 

WHILE AS = "" 




’ A$ = INKEYS 
’ WEND 
* CLS 
RETURN 

getstartadd: ’read 18 bit stop address from 3 ports and convert to number 
’only 16 bit stop address from 2 ports for 64 k sample memory 
stopadd = 0 
’bits 16 and 17 

REM saddl = INP(physicaladdress(i) + 2) AND &H3 
REM saddl = saddl * 64 * 1024 
REM stopadd = saddl 
’bits 8 thru 15 

saddl = lNP(physicaladdress(r) + 1) 
saddl = saddl * 256 
stopadd = stopadd + sadd 1 
’bits 0 thru 7 

stopadd = INP(physicaladdress(I)) + stopadd 
’Calculate starting address by subtracting samples 
starts(I) = stopadd - SAMPLES - 4 
’Starting address is never negative, so fix it if needed. 

IF starts(I) < 0 THEN stans(I) = starts(I) + 65536 
RETURN 

’GETDAT and gdl both get a single value from the recorder memory, 
’this value is returned in the variable DAT. DAT is a 
’number that ranges from -2048 to +2047. 

’GETDAT will initialize die mur.cy pointer with the value in BRDST. 
’However, this memory pointer need only be initialized one time per 
’record transfer. Then the calling routine can continuously call gdl 
’and the recorder memory pointer will auto increment after each call. 

getdat: 

IF brdst > memmax THEN brdst = brdst - memmax 
IF brdst < 0 THEN brdst = brdst + memmax 
REM str = INT(brdst / 65536) 

REM brdst = brdst - (str * 64 * 1024) 
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strhi = INT(brdst / 256) 
strlo = brdst - (strhi * 256) 

OUT physicaladdress(I), strlo 
OUT physicaladdress(I) + 1, strhi 
REM OUT physicaladdress(i) + 2, str 

gdl: 

A = INP(physicaladdress(I) + 3) 

B = INP(ph)'sicaladdress( I ) + 4) 

B = (B AND &HF) * 256 
DAT = (A + B - 2048) * -1 
RETURN 


PLOT: ’This rcaiinc will plot the arrays to the screen. 

SC KLIN 2 ’2 for cga, 3 for ties and 

’9 for ega. (Only requires for 
’plotting) 

WINDOW (0, 2048)-(SAMPLES, -2048) 

’plot data on screen 
FOR I = 1 TO 2 


GOSUB getstanadd ’Determines the stanmg point of 
’the record in memory 

brdst = atarts(l; ’Initializes memory pointer and gets 
GOSUB getdat 'first data point 

PSET 0), DaT) 

FOR , = 1 TO SAMPLES 

GOSUB gdl ’gets sample in dal and auto inc. ptr 
LINE -(J, DAT) 

NEXT J 
NEXT I 
RETURN 


LOAD. T^tad data into an array 


FOR I = 1 TO 2 

GOSUB gctstaiadd 


’Initializes boards for data 




brdst = starts® 'downloading and gets first 
GOSUB gctdat 'sample 

FOR J = 1 TO (2 * SAMPLES) STEP 2 

GOSUB gd 1 'Gets sample and increments pointer 

IF I = 1 THEN 

VRAW(J) = DAT 
VRAW(J + 1) = 0! 

END IF 
IF I = 2 THEN 

PRAW(J) = DAT 
PRAW(J + 1) = 0! 

END IF 


NEXT J 
NEXT I 
RETURN 


HEADER: 

'This subroutine collects run information for a file header and 
'data processing 

CLS 'Clear screen 

LOCATE 2, 10 'Locates text 

PRINT "RUN DATA INPUT" 

LOCATE 4, 20 

PRINT "1. Input a Run Number consisting of the Model" 
LOCATE 5, 25 

PRINT "Number and the Angle of Attack in the form XYY." 
INPUT RUNNOS 
LOCATE 6, 25 

PRINT "RUN NUMBER - "; RUNNOS 
LOCATE X, 20 
PRINT "2 Input Date' 

INPUT DA ITS 
LOCATE d, 25 
PRINT "DATE- . DA PIS 
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LOCATE 11, 20 
PRINT ''3. Input Time" 

INPUT TYMES 
LOCATE 12, 25 
PRINT "TIME-"; TYMES 
LOCATE 14, 20 

PRINT "4. Input Model Description" 

INPUT MODELS 
LOCATE 15, 25 
PRINT "MODEL-"; MODELS 
LOCATE 17, 20 

PRINT "5. Inpui Angle of Attack (+/-) in Degrees" 

INPUT AO A 
LOCATE 18, 25 

PRINT "ANGLE OF ATTACK-"; AOA; "DEGREES" 

LOCATE 20, 20 

PRINT “6. Input Mean Dynamic Pressure in CM of Water" 

INPUT PDYN 
LOCATE 21, 25 

PRINT "MEAN DYNAMIC PRESSURE-”; PDYN; "CM-WATER" 
LOCATE 23, 20 

PRINT "7. Input Indicated Oscillation Frequency (Hz)" 

INPUT OMEGA 
LOCATE 24, 25 
CLS 

LOCATE 3, 20 

PRINT OSCILLATION FREQUENCY-"; OMEGA; "Hz" 
LOCATE 5, 20 

PRINT "8. Input Barometric Pressure (In He)" 

INPUT PREF 
LOCATE 7, 25 

PRINT "BAROMETRIC (REF) PRESSURE-"; PREF; "IN HG" 
LOCATE 9, 20 

PRINT "9. INPUT Frecstream Temperature (F)" 

INPUT TEMP 
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LOCATE 10, 25 

PRINT "FREESTREAM TEMPERATURE-"; TEMP; "F" 

RETURN 


HOME: 

’This subroutine advances the SCANIVALVE to HOME 
FOR I = 1 TO 50000 
OUT &H2A8, 2 
NEXT I 

OUT &H2A8, 0 

RETURN 


PORTNUM: 

’This file contains the text of the subprogram PORTNUM in TKDATA 

’This subprogram reads the BCD output of the SCANIVALVE, 
’(PORTNO) 

PORTNO = 255 - INP(&H2A9) 

RETURN 


PORT: 

’This file contains the text of the subprogram PORT in TKDATA 
’This subprogram advances the SCANIVALVE one port 
FOR I = 1 TO 1000 
OUT &H2A„, 1 
NEXT 1 

OUT &H2A8, 0 

RETURN 


MEAN: ’Determines the mean values of the two arrays PRAW and VRAW, PBAR 
’and VBAR, and replaces the arrays with the mean values subtracted 
VBAR = 0! 

PBAR = 0! 

FOR 1 = 1 TO (2 * SAMPLES) STEP 2 

VBAR = CDBL(VRAW(I) * (VSCALE / 2(U8) + VBAR) 

PBAR = CDBL(PRAWO) * (PSCALE / 2(S48) + PBAR) 
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NEXT I 

VBAR = VBAR / SAMPLES 

PBAR = PBAR / SAMPLES 

FOR I = 1 TO (2 * SAMPLES) STEP 2 

VRAW(I) = VRAW(I) * (VSCALE / 2048) - VBAR 
PRAW(I) = PRAW(I) * (PSCALE / 2048) - PBAR 
NEXT I 
RETURN 

FFTPRAW: ’Replaces PR AW with the FFT of PR AW, (in pseudo-complex form) 
n = 2 * SAMPLES 
J = 1 

FOR I = 1 TO n STEP 2 
IF J > I THEN 

TEMPR = PRAW(J) 

TEMPI = PRAW(J + 1) 

PRAW(J) = PRAW(I) 

PRAW(J + 1) = PRAW(I + 1) 

PRAW(I) = TEMPR 
PRAW(I + 1) = TEMPI 
END IF 
m = n / 2 

1 IF m > 2! AND J > m THEN 

J = J - m 
m = m / 2 
GOTO 1 
END IF 

J = J + m 

NEXT I 

mntax = 2 

2 IF n > mrnax THEN 

ISTEP = 2 * mm ax 
THETA = 6.283185307# / mmax 
WPR = -2! * SIN(.5 * THETA) A 2 
WPI = SIN(THETA) 








WR = 1! 

WI = 0! 

FOR m = 1 TO mmax STEP 2 

FOR I = m TO n STEP ISTEP 
J = I + mmax 

TEMPR = WR * PRAW(J) - WI * PRAW(J + 1) 
TEMPI = WR * PRAW(J + 1) + WI * PRAW(J) 
PRAW(J) = PRAW(I) - TEMPR 
PRAW(J + 1) = PRAW(I + 1) - TEMPI 
PRAW(I) = PRAW(I) + TEMPR 
PRAW(I + 1) = PRAW(I + 1) + TEMPI 
NEXT I 

WTEMP = WR 

WR = WR * WPR - WI * WPI + WR 
WI = WI * WPR + WTEMP * WPI + WI 
NEXT m 

mmax = ISTEP 
GOTO 2 

END IF 

RETURN 


XFERINPUT: This subroutine loads transfer function values in the standard 
’format from an external f.'.e XFER.FCN (which MUST be in the 
’same director) as TKDATA). 

OPEN "XFER.FCN" FOR INPUT AS #1 
FOR I = 1 TO 2048 

INPUT #1, XFERJ, 

NEXT I 
RETURN 


IFFTRAW; ’Replaces PRAW with die IFFT of PRAW, (in pseudo-complex form) 
n = 2 * SAMPLES 
J = 1 

FOR I = 1 TO n STEP 2 
IF J > I THEN 
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TEMPR = PRAW(J) 

TEMPI = PRAW(J + 1) 

PRAW(J) = PRAW(I) 

PRAW(J + 1) = PRAW(I + 1) 

PRAW(I) = TEMPR 
PRAW(I + 1) = TEMPI 
END IF 
m = n / 2 

3 IF m > 2! AND J > m THEN 

J = J - m 
m = m / 2 
GOTO 3 
END IF 

J = J + m 

NEXT I 

mmax = 2 

4 IF n > mmax THEN 

ISTEP = 2 * mmax 
THETA = -6.283185307# / mmax 
WPR — -2! * SIN(.5 * THETA) A 2 
WPI = SIN(THETA) 

WR = 1! 

WI = 0! 

FOR m = 1 TO mmax STEP 2 

FOR I = m TO n STEP ISTEP 
J = I + mmax 

TEMPR = WR * PRAW(J) - W’l * PRAW(J + 1) 
TEMPI - WR * PRAW(J + 1) + WI * PRAW(J) 
PRAW r (J) = PRAWd) - TEMPR 
PRAW(J + 1) = PRAW(I + 1) - TEMPI 
PRAW(I) = PRAW(I; + TEMPR 
PRAW(1 tl) = PRAWlI + 1; + TEMPI 
NEXT I 

W'TEMP = WR 

WR = WR * WPR - WI * WPI + W ; R 

X 1 




WI = WI * WPR + WTEMP * WPI + WI 
NEXT m 

mmax = ISTEP 
GOTO 4 

END IF 

FOR 1 = 1 TO n 

PRAW(I) = PRAW(I) / SAMPLES 
NEXT I 
RETURN 


RIRTHETA: ’This subroutine converts an array from real-imaginary pairs to 
’r-thcta pairs. 

FOR 1 = 1 TO 2047 STEP 2 
IF I = 1 THEN 
R = 0! 

THETA = 0! 

ELSE 

R = SQR(PRAW(I) a 2 + PRAW(I + 1) A 2) 

IF PRAW(I) = 0! THEN 

THETA = 1.570796327# * SGN(PRAW(I + 1)) 

ELSE 

THETA = ATN(PRAW(I + 1) / PRAWfJ); 

END IF 
END IF 

PRAW(I) = R 
PRAW(I +U = THETA 
NEXT I 
RETURN 


RTHETARI: ’This subroutine converts an array from r-thcta pairs to rcal- 
’imaginary pairs. 

FOR m = 1 TO 2047 STEP 2 

R = PRAW(m) * COSCPRAWtm + 1)) 

I = PRAWfm) * SINfPRAWfm + 1)) 
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PRAW(m) = R 
PRAW(m + 1) = I 
NEXT m 


RETURN 


DIVIDE: ’This subroutine replaces PRAW with PRAW/XFER. 
FOR m = 1 TO 2047 STEP 2 

PRAW(m) = PRAW(m) / XFER(m) 

PRAW(m + 1) = PRAW(m + 1) + XFER(m + 1) 
NEXT m 
RETURN 


ARRAY: ’This subroutine loads PBAR and VBAR into vectors PMEAN and VMEAN, 
’and PRAW and VRAW into arrays P and V. (It truncates the 1024 
’records in PRAW and VRAW into 1000 records, performs an ensemble 
’average over the 10 ensembles and loads a 100 sample array ) 

PMEAN(PORTNO) = PBAR 
VMEAN(PORTNO) = VBAR 
FOR k = 1 TO 100 
PEA = 0! 

VEA = 0! 

FOR J = 0 TO 9 

PEA = PRAW((2 * k) - 1 + (J * 100)) 

VEA = VRAW((2 * k) - 1 + (J * 10(0) 

NEXT J 

P(PORTNO, k) = PEA 
V(PORTNO, k) = VEA 
NEXT k 
RETURN 


NEGALPHA: This subroutine halts data acquisition to allow the angle of 
’attack to be changed from positive alpha to negative alpha, 

CLS 
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BEEP 

LOCATE 8, 10 

PRINT ”1.) Program execution has been interrupted to allow angle of' 

LOCATE 10, 15 
PRINT “attack change." 

LOCATE 15, 10 

PRINT "3.) Change ALPHA to MINUS ALPHA then press any key to continue” . 

LOCATE 17, 15 
PRINT "data collection." 

AS = "" 

WHILE AS = "" 

AS = INKEYS 

WEND 

CLS 

LOCATE 10, 30 

PRINT "Data collection continuing." 

RETURN 


DISK: ’This subroutine creates a file AtRUNXXX.DAT, (where XXX is the 

’RUN NUMBER), and writes all HEADER DATA, the mean pressure and 
’velocity vectors, PMEAN and VMEAM and the pressure and velocity 
’arrays, P and V to that file. 

CLS 

LOCATE 10, 15 

PRINT "INSERT A FORMATTED DISK IN DRIVE A:" 

LOCATE 15, 15 

PRINT "PRESS ANY KEY TO WRITE RUN DATA TO FILE" 

AS = "" 

WHILE AS = "" 

AS = INKEYS 

WEND 

CLS 

LOCATE 10, 10 
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PRINT "Run data is being written to file A:RUN" + RUNNOS + ".DAT" 
’Open the file: 

OPEN "A:RUN" + RUNNOS + ".DAT" FOR OUTPUT AS #2 
’Write header data: 

WRITE #2, RUNNOS, DAITES, TYMES, MODELS: ’Strings! 
WRITE #2, AOA: ’degrees 
WRITE #2, PDYN * 2.048: ’psf 
WRITE #2, OMEGA: ’hertz 
WRITE #2, PREF * 70.73: ’psf 
WRITE #2, TEMP + 459.7: ’deg Rankine 
’Next write pressure and velocity data: 

FOR I = 1 TO 46 

WRITE #2, PMEAN(I): ’Mean pressure vector 
NEXT I 

FOR I = 1 TO 46 

WRITE #2, VMEAN(I): ’Mean velocity vector 
NEXT I 

FOR I = 1 TO 46 

FOR J = 1 TO 100 

WRITE #2, P(I, J): ’Pressure array 
NEXT J 
NEXT I 

FOR I = 1 TO 46 

FOR J = 1 TO 100 

WRITE #2, V(I, J): ’Velocity array 
NEXT J 
NEXT I 
CLS 
BEEP 

LOCATE 10, 10 

PRINT " RUN COMPLETE AND DATA WRITTEN TO FILE" 


RETURN 





APPENDIX B 


The following is the data reduction program, REDUCE.FOR, which is written in 
FORTRAN. 


PROGRAM REDUCE 
c 

c ******** ********* ********* ********* ********* ********* ********* ** 

c 

c DATA REDUCTION PROGRAM FOR THESIS. WRITTEN BY 
c LT G. A. CHAMBERLAIN. JR., USN, JANUARY 1990. 
c CURRENT UPDATED VERSION, 11 JULY 1990 
c 

c THIS VERSION IS INTENDED FOR USE ON THE PC IN BLDG 230. 

C PROGRAM "REDUCE.FOR’ SHOULD BE INSTALLED IN THE ROOT DIRECTORY OF 
c DRIVE C:, AND FILE "XOVRC DAT" SHOULD BE IN THE ROOT DIRECTORY OF 
c DRIVE D:. THIS PROGRAM WILL ASK FOR THE APPROPRIATE RUN NUMBER, 
c ONCE INPUT, THE PROGRAM WILL READ FILE "RUNXXX.DAT" FROM DRIVE 
c C:\DATA AND CUE THE USER AS TO THE PROGRESS OF THE DATA REDUCTION, 
c OUTPUT FILES ARE PUPXXX.DAT, PLOXXX.DAT, VUPXXX.DAT, VLOXXX.DAT, 
c CPUXXX.DAT AND CPLXXX.DAT. ALL OF THESE WILL BE WRITTEN TO DRIVE C: 
c IN DIRECTORY C:\DATA\filciiaiTie.exi. ADDITIONALLY, FILES CL.DAT, 
c CD.DAT AND CMC4.DAT (WHICH MUST EXIST ALREADY IN C:\DATA\) WILL 
c BE APPENDED WITH THE CL, CD AND CMC4 INFORMATION FOR THE AOA 
c USED IN THAT EXPERIMENTAL RUN. 
c 
c 

c FURTHER DATA MANIPULATION WITH SURFER AND GRAPHER IS FACILITATED 
c BY SETTING THE DEFAULT DRIVES IN THE "ENVIRONMENT" MENUS OF THESE 
c PROGRAMS TO LOOK FOR INPUT AND WRITE OUTPUT TO C:\DATAV ALSO, 


86 










c DATA FROM EACH RUN CAN BE EASILY COLLECTED AND COPIED TO THE 
c APPROPRIATE DISK FOR STORAGE BY USING THE DOS COPY COMMAND, AS 
C FOLLOWS: " \> COPY ???123 * A: IT IS IMPORTANT TO NOTE THAT 
c AS LONG AS ALL FILES THAT ARE CREATED HAVE THE SAME NAME STRUCTURE, 
- c ie, ABC123.exe, THEN A SINGLE COPY COMMAND WILL CAUSE THE COMPUTER 

c TO RETRIEVE ALL FILES THAT ARE SO NAMED. 

« c 

c *** KEY TO LOGICAL UNIT NUMBERS *** 
c 7...RUNXYY.DAT, INPUT DATA FILE 

c 8-.XOVRC.DAT, CONTAINS PRESSURE PORT LOCATION DATA 
c 9...CPUXYY.DAT, OUTPUT FILE CONTAINING CPupper vs TIME AND LOCATION 
c 10..CL.DAT, OUTPUT FILE CONTAINING CL vs TIME AND LOCATION 
c 11 ...CD.DAT, OUTPUT FILE CONTAINING CD 
c 12...CMC4.DAT, OUTPUT FILE CONTAINING CMC4 

c 19...CPLXYY.DAT, OUTPUT FILE CONTAINING CPlower vs TIME AND LOCATION 
c 20...PUPXYY.DAT, OUTPUT FILE W/ UPPER SURFACE PRESSURE PROFILE 
c 21 ...PLOXYY.DAT, OUTPUT FILE W/ LOWER SURFACE PRESSURE PROFILE 
c 22...VUPXYY.DAT, OUTPUT FILE W/ UPPER SURFACE VELOCITY PROFILE 
c 23...VLXYY.DAT, OUTPUT FILE W/ LOWER SURFACE VELOCITY PROFILE 
c 

a £******** ********* ********* ********* ********* ********* ********* ** 

c 

IMPLICIT DOUBLE PRECISION (A-H,L,0-Z; 

DIMENSION PRM(46j.VM(46j,CPUPR(23,100),CPLWR(23,100j,X(23;, 

\ CPU2(23),CPL2(23),PRES(46,100), 

\ VELOC(46,100), V2D(100), CL(100), CD(100), CMC4(100; 

CHARACTER XXX*3, RUNID*10 

WRITE (6,*j ’RUN IDENTIFICATION NUMBER CONSISTS OF THREE’ 

WRITE (6,*) ’DIGITS, XYY, WHERE X = MODEL SHAPE (1,2 or 3)' 

WRITE (6,*; ’AND YY = ANGLE OF ATTACK IN DEGREES.’ 

DO 321 I = 1,K) 

WRITE (6,*) 

. 321 CONTINUE 

4 WRITE (6,*) ’INPUT RUN IDENTIFICATION NUMBER XYY ?’ 

READ (5,1000) XXX 
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1000 FORMAT (A3) 

DO 999 I = 1,30 
WRITE (6,*) 

999 CONTINUE 

RUNID = ’RUN’//XXX//’.DAT’ 

5 WRITE (6,*) ’ENSURE DATA DISK CONTAINING FILE \RUNID,’ IS ON’ 

WRITE (6,*) ’HARD DISK AT C:\DATA\ . REDUCE FOR WILL LOOK FOR’ 
WRITE (6,*) RUNID,’ AT THAT LOCATION. IF NOT, PRESS CTRL-C NOW’ 
WRITE (6,*) ’AND COPY ’.RUNID,’ TO C:\DATA\ THEN BEGIN AGAIN.’ 

WRITE (6,*) 

WRITE (6,*) 

WRITE (6,*) ’OUTPUT FILES WILL BE PUPXYY.DAT, PLOXYY.DAT, ’ 

WRITE (6,*) ’VUPXYY.DAT, VLOXYY.DAT, CPUXYY.DAT AND CPLXYY.DAT.’ 
WRITE (6,*) ’ALL OF THESE WILL BE WRITTEN TO C:NDATA\ AND MAY’ 
WRITE (6.*) ’BE RETRIEVED TO THE "A" DISK BY THE DOS COMMAND’ 
WRITE (6,*) ’... COPY C :\D AT AY??? 123.* A: ... ’ 

WRITE (6,*) ’CL, CD AND CMC4 FILES WILL ALSO BE AMENDED TO’ 

WRITE (6,*) ’THIS NEW AOA.’ 

DO 6 I = 1,6 
WRITE (6*) 

6 CONTINUE 

PAUSE ’STRIKE <ENTER> TO CONTINUE’ 

WRITE (6,*) ’PROCEEDING...’ 

OPEN (UNIT = 7, FILE = ’C:\DATA\RUN’//XXX//’.DAT’, 

/ STATUS = ’OLD’) 

OPEN (UNIT = 8, FILE = ’D:\XOVRC.DAT’, STATUS = ’OLD’) 

OPEN (UNIT = 9, FILE = ’C:\DATA\CPU7/XXX//'.DAT’, 

/ STATUS = ’NEW’) 

OPEN (UNIT = 19, FILE = ’C:\DATA\CPL’//XXX//’.DAT’, 

/ STATUS = ’NEW’) 

OPEN (UNIT = 20, FILE = ’C:\DATA\PUP’//XXX//’.DAT’, 

/ STATUS = ’NEW’) 

OPEN (UNIT = 21, FILE = ’C:\DATA\PLO’//XXX//’.DAT’, 

/ STATUS = ’NEW’) 

OPEN (UNIT = 22. FILE = ’C:\DATA\VUP’//XXX//’.DAT’. 
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/ STATUS = ’NEW’) 

OPEN (UNIT = 23, FILE = ’C:\DATAWLOV/XXX//’.DAT’, 

/ STATUS = ’NEW’) 

READ (7,*) 

READ (7,*) A 
40 CONTINUE 

PI = 3.141592653589 
ALFA = A * PI / 180.0 
READ (7,*) Q 
READ (7,*) 

READ (7,*) P 
READ (7,*) T 

WRITE (6,*) ’READING MEAN PRESSURES’ 

DO 60 I = 1,46 
READ (7,*) PRM(I) 

60 CONTINUE 

WRITE (6,*) ’READING MEAN VELOCITIES’ 

DO 70 I = 1,46 
READ (1*) VM(I) 

70 CONTINUE 

WRITE (6,*) ’READING DYNAMIC PRESSURE CORRECTIONS & COMPUTING’ 
WRITE (6,*) ’PRESSURE DISTRIBUTION AND COEFFICIENTS.’ 

DO 80 I = 1,23 
DO 90 J = 1,100 
READ (7,*) PRD 
PRES (I ,]) = PRD + PRM(I) 

CPUPR(U) = l+(PRES(U)/0) 

90 CONTINUE 
80 CONTINUE 
DO 10(1 I = 1.23 
K = I + 23 
DO 1! 0 J = 1.100 
READ (7*) PRD 
PRES (K.J) = PRD + PRM (K) 

CPLWR(I.J) = 1 + (PRESIKJ) IQ) 
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110 CONTINUE 
100 CONTINUE 

WRITE (6,*) ’READING DYNAMIC VELOCITY CORRECTION & COMPUTING’ 
WRITE (6,*) ’VELOCITY DISTRIBUTION OVER THE AIRFOIL.’ 

IX) 120 I = 1,46 
DO 130 J = 1,100 
READ (7,*) VD 
VELOC(LJ) = VD + VM(I) 

130 CONTINUE 
120 CONTINUE 
DO 140 I = 1,23 
READ (8,*) X(I) 

140 CONTINUE 
DO 998 I = 1,30 
WRITE (6,*) 

998 CONTINUE 

WRITE (6,*) ’WRITING PRESSUREA'ELOCITY ARRAYS TO HARD DISK:’ 
WRITE (6,*) ’UPPER SIDE PRESSURE...’ 

DO 997 I = 1,23 
DO 996 J = 1,100 

WRITE (20,6996) X(I),J/100.0,(P+PRES(I,J)) 

996 CONTINUE 

997 CONTINUE 
CLOSE (20) 

WRITE (6,*)’ ...LOWER SIDE PRESSURE...’ 

DO 995 I = 1,23 
K = I + 23 
DO 994 J = 1,100 

WRITE (21,6996) X(I),J/100.0,(P+PRES(K,J» 

994 CONTINUE 

995 CONTINUE 
CLOSE (21) 

WRITE (6,*) ’UPPER SIDE VELOCITY...’ 

DO 222 I = 1,23 
DO 223 J = 1,100 
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WRITE (22,6996) 1/23.0, J/100.0, VELOC(l,J) 

223 CONTINUE 
222 CONTINUE 

CLOSE (22) 

WRITE (6,*) ’ ...LOWER SIDE VELOCITY...’ 

DO 224 I = 1,23 
K = I + 23 
DO 225 I = 1,100 

WRITE (23,6996) U23.0, J/100.0, VELOC(K,J) 

225 CONTINUE 

224 CONTINUE 
CLOSE (23) 

6996 FORMAT (1X,F20.5,1X,F20.5,5X,F20.5) 

WRITE (6 *) 'WRITING UPPER SURFACE CP FILE TO HARD DISK:’ 

WRITE (6,*) ’(DONT FORGET THAT THIS WILL BE NEGATIVE CP)’ 

DO 150 I = 1,23 
DO 155 J = 1,100 

WRITE (9,6000) X(I),J/100.0,(-CPUPR(U)) 

155 CONTINUE 
150 CONTINUE 

WRITE (6,*) 'WRITING LOWER SURFACE CP FILE TO HARD DISK:’ 

WRITE (6,*) ’(ALSO NEGATIVE CP)’ 

DO 160 I = 1,23 
DO 165 J = 1,100 

WRITE (19,6000J X(I),J/100.0,(-CPLWR(I,J)) 

165 CONTINUE 
160 CONTINUE 

6000 FORMAT (IX, F6.3.F20.5.E15.5) 

CLOSE (9) 

CLOSE (19) 

WRITE (6,*) ’PRESS CONTROL-C NOW TO AVOID SPLINE INTERPOLATION 
WRITE (6,*) ’PRESSURE, VELOCITY AND Cp FILES ARE ALREADY SAVED. 
WRITE (6,*) 

WRITE (6,*) ’BEGINNING PERFORMANCE COEFFICIENT CALCULATIONS’ 
DO 777 J = 1,100 
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CMC4(J)=0.0 
CL(J)=0.0 
DO 778 I = 1,22 


DCL=(CPLWR(1,J)-CPUPR(1,J))* ,5*(X(I+ l)-X(l-l)) 
CL(J)=CL(J)+DCL 

CMC4(J)=CMC4(J)+<0.25-X(I))*DCL 
778 CONTINUE 

CLF=(CPLWR(U)-CPUPR(1,J))*.0125 

CMF=CLF*(0.25-X(1)) 

CLB=(CPLWR(23,J)-CPUPR(23,J))*.075 
CMB=CLB*(0.25 X(23)) 

CL(J)=(CL(J)+CLF+CLB)*COS(ALFA) 

CD(J)=CL(J)*TAN(ALFA) 

CMC4(J)=CMC4(J)+CMF+CMB 
111 CONTINUE 

c ****** STEADY FLOW LOOP, PRINTS AVG CL, CD, CM ******** 

C ****** VALUE TO SCREEN FOR IMMEDIATE FEEDBACK ********* 

c AVG1 = 0.0 

c AVG2 = 0.0 

c AVG3 = 0.0 

c DO 1234 I = 1,100 

c AVG1 = AVG1 + CL(I) 

c AVG2 = AVG2 + CD(1) 

c AVG3 = AVG3 + CMC4(I) 

c WRITE (6,8888) A,I,CL(I),CD(I),CMC4(I) 

cl 234 CONTINUE 

c A VG1 = AVG 1/100.0 

c AVG2 = AVG2/100.0 

c AVG 3 = AVG3/100.0 

c WRITE (6,8889) A, AVG1, AVG2, AVG3 

c8888 FORMAT (IX,F5.1,14,3F10.4,/) 

c8889 FORMAT (IX,’AVG VALUES AT \F8.1,’ DEGREES AOA:’,/,lX,’CL 
c / FI0.4,’ CD = \F10.4,' CMC4 = ’.F10.4,/) 

C **»*****END OF STEADY FLOW LOOP *********************** 

OPEN (UNIT = 10, FILE = ’C:\DATA\CL.DAT’, STATUS = ’OLD’) 
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DO 10 I = 1,4000 

READ (10,* ,ERR= 1 ,IOSTAT = J10) Z 
IF (J10.LT.0) GO TO 1 
10 CONTINUE 

1 BACKSPACE (10) 

DO 220 J = 1,100 

WRITE (10,6100) A, J, CL(J) 

220 CONTINUE 

OPEN (UNIT = 11, FILE = ’C:\DATA\CD.DAT’, STATUS = ’OLD’) 

DO 20 I = 1.4000 

READ (11 *,ERR=2,IOSTAT = Jll) Z 
IF (J11.LT.0) GO TO 2 
20 CONTINUE 

2 BACKSPACE (11) 

DO 230 J= 1,100 

WRITE (11,6100) A, J, CD(J) 

230 CONTINUE 

OPEN (UNIT = 12, FILE = ’C:\DATA\CMC4.DAT’,STATUS = ’OLD’) 

DO 30 I = 1,4000 

READ (12,*,ERR=3,IOSTAT = J12) Z 
IF (J12.LT.0) GO TO 3 
30 CONTINUE 

3 BACKSPACE (12) 

DO 240 J = 1,100 

WRITE (12,6100) A, J, CMC4(J) 

240 CONTINUE 

6100 FORMAT (lx, F6.2, I4,E15.5) 

CLOSE (7) 

CLOSE (8) 

CLOSE (10) 

CLOSE (11) 

CLOSE (12) 

WRITE (6,*) ’DATA REDUCTION COMPLETE.’ 

£******** ********* ********* ********* ********* ********* ********* ** 
END 
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